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SUMMARY 


An analysis was made which permitted calculation of the effective- 
ness of oxide coatings in retarding the transient flow of heat into tur- 
bine rotor blades when the combustion-gas temperature of a turbojet 
engine is suddenly changed. In order to use the analysis, it was nec- 
essary to know the coating thickness, the coating thermal conductivity, 
and the turbine flow conditions. An experimental investigation was made 
to determine the thickness of coating that will build up on the turbine 
blades and to check the reliability of the analysis for calculating the 
turbine-blade temperature response when the engine operating conditions 
are suddenly changed. Two fuel additives, SF-99 silicone oil and 
tributyl borate, were mixed with the MIL-F-5624A (grade JP-4) fuel to 
form coatings of silicon dioxide and boric oxide, respectively, on the 
turbine blades. The experimental investigation consisted in recording 
the transient temperature at several locations on the turbine blade while 
the engine was both accelerated and decelerated between the speeds of 
8000 and 11,500 rpm. The engine was operated with straight JP-4 fuel as 
well as JP-4 fuel mixed separately with each of the fuel additives. The 
uncoated-blade temperature response wes compared with coated-blade 
temperature responses to determine the effectiveness of the two coatings. 


The very thin silicon dioxide and boric oxide coatings that formed 
on the turbine rotor blades (approximately 0-001 in.) resulted in a 
negligible effect on the blade temperature response. Calculations using 
the analysis of this report showed that for this coating thickness a 
coating conductivity of 0.008 Btu per hour per foot per ОР would be 
needed to produce a lag of only 100° F in blade temperature relative to 
the uncoated blade 60 seconds after the engine operating conditions were 
changed from 8000 to 11,500 rpm. Мо oxide coatings are currently avail- 
able with thermal conductivities this low. А comparison of experimental 
with calculated transient uncoated-turbine-blade temperature showed that 
the blade-metal temperature calculated from the analysis of this report 
was а maximum of 409 F lower than the experimental turbine rotor-blade 
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temperature during an engine acceleration from 8000 to 11,500 rpm. The 
use of fuel additives had no measurable effect on the engine thrust. 


INTRODUCTION 


The combustion-gas temperature of & turbojet engine is limited to 
а value which will result in a reasonable turbine-blade life. If it 
were possible to increase the combustion-gas temperature for a short 
period of time (60 sec) without the turbine-blade temperatures exceeding 
& limiting value, the thrust of a turbojet engine could be increased for 
short periods of time, provided the compressor does not surge as the 
engine speed and temperature are increased. A low thermal conductivity 
coating on the outside of the turbine blades would retard the flow of: 
heat into the blades and thereby lengthen the time during which the 


extra thrust is available. Reference 1 shows that the steady-state heat- 


flow rate into water-cooled rocket motor walls was decreased from 30 to 
50 percent by 1/8-inch-thick coatings that were applied dynamically by 
condensing the coating material from the products of combustion on the 
rocket motor walls. The analytical study reported in reference 2 shows 
that a 0.015-inch coating having a thermal conductivity of 0.25 Btu per 
hour per foot per OF would be effective in appreciably reducing the 
trailing-edge temperature of water-cooled turbine blades. Strong cool- 
ing existed in both references 1 and 2. 


The results of an experimental investigation made to determine the 


effect of dynamically applied coatings of silicon dioxide and boric oxide 


on the transient temperature of a single, stationary, uncooled turbine 
blade made of S-816 alloy are reported in reference 3. Іп this refer- 
ence transient blade temperature data were obtained on both the coated 
and the uncoated blade when the gas temperature was suddenly increased. 
The results showed a negligible effect of the coating on the rate of — 
blade temperature response. 


Although the coatings used in reference 3 were ineffective in 
decreasing the temperature-response rate of a stationary turbine blade, 
it is still desirable to be able to calculate the coating properties 
required for effectively retarding the transient flow of heat into gas- 
turbine blades when engine operating conditions are changed. Before 
these coating properties can be calculated, it is necessary to determine 
both the accuracy of the calculation method and the coating thicknesses 
that can be expected to form on the turbine blades. In order to deter- 
mine these factors, both analytical and experimental investigations were 
made at the NACA Lewis laboratory. 


The objectives of the analytical investigation were to obtain (1) a 
method of evaluating experimental data to determine the effectiveness of 
different coating materials in decreasing the temperature-response rate 
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of turbine blades and (2) an analytical method for determining the time- 
temperature response of a coated turbine blade. The objectives of the 
experimental investigation were to determine (1) the coating thickness 
that will occur on the turbine blades іп a full-scale turbojet engine 
when the coating material is either silicon dioxide or boric oxide, (2) 
the effect of these coatings on the transient turbine-biade temperatures, 
(3) the accuracy with which the transient turbine-blade temperatures can 
be calculated from the analysis presented in this report, and (4) the 
effect of the coating on the engine parts and on the over-all engine 
performance. 


The experimental investigation was made on & production turbojet 
engine the turbine blades of which were instrumented with thermocouples 
in the leading-edge, midchord, and trailing-edge regions. Transient 
blade temperature data for both engine acceleration and deceleration 
were obtained for the engine speed range from 8000 to 11,500 rpm. These 
speeds were chosen because they represented the extremes of combustion- 
gas temperature that could be obtained in the operation of the engine 
used in this investigation. 


Two fuel additives, SF-99 silicone oil and tributyl borate, were 
used in the investigation to form coatings of silicon dioxide and boric 
oxide, respectively. Transient turbine-blade temperature data were 
obtained first with straight MIL-F-5624A (grade JP-4) fuel to obtain 
uncoated-turbine-blade data, and then similar data were obtained when 
each of the fuel additives was mixed with the fuel. This time- 
temperature data permitted the calculation of time constants (the time 
required for the blade temperature to increase or to decrease approxi- 
mately 63 percent of the over-all difference between the initial and 
final blade temperatures), which were used to evaluate the effectiveness 
of the coatings on the blade temperature response. Data were taken dur- 
ing both engine acceleration and deceleration to provide additional 
transient temperature data for comparison with the analysis of this 
report. 


After each fuel additive was used, the engine parts which were 
exposed to the combustion gases were inspected to determine the effect 
each coating might have on the engine parts. The effect of each coating 
on the engine performance was also observed during the investigation. 


ANALYSIS 


The strength of most turbine-blade metals decreases very rapidly 
with increases in the blade-metal temperature. Therefore, blade temper- 
ature increases even of short duration can cause blade failure. For 
this reason, if the gas temperature is to be raised above the design 
value, it would be convenient to be able to predict the rate of increase 
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of the blade temperature under the influence of the increased gas tem- 
perature. Then the time for the blade temperature to reach a critical 
value could be calculated. Also, if insulating coatings are used on 
the blade surface to retard the flow of heat into the blade, a method 
should be available to predict the effectiveness of these coatings in 
decreasing the rate of blade-metal temperature response to changes in 
gas temperature. Such an analysis of transient blade-metal temperatures 
at the midchord of both coated and uncoated blades is made in this sec- 
tion. A list of the symbols used in this and other sections of the 
report is presented in appendix A, and & similar analysis for the tran- 
sient metal temperatures at the leading or trailing edges of turbine 
blades is given in appendix B. 


Uncoated Blade `` 


For the analysis of the transient blade-metal temperatures at the 
midchord location, a blade section was chosen at a sufficient distance 
from the blade base, so that there would be no spanwise heat flow into 
the section. An element of this section was chosen at the midchord 
Location where there would be no chordwise heat flow across the element 
boundaries. Such an element of a blade section is shown in figure 1. 
The length of the element is designated a (the distance between the 
pressure and the suction surface at the chosen location). The width of 
the element is b. The opposite sides of the element are assumed 


parallel. 


The analysis considers two cases: one in which the blade-metal 
thermal conductivity is infinite, and one in which the blade-metal 
conductivity is finite. Considering the blade thermal conductivity 
infinite greatly simplifies the analysis; but since the turbine-blade- 
metal thermal conductivity is low (about 10 to 15 Btu/(ft)(hr)(9F)), the 
effect of the blade thermal conductivity may be significant. For this 
reason, both cases are considered in this analysis, and the blade tem- 
peratures calculated by both theoretical equations are compared to 
determine the error introduced by considering the metal thermal con- 
ductivity infinite. 


Case 1 - infinite blade-metal thermal conductivity. - For this case 


there will be no temperature gradients across the blade section. Con- 
sider the strip of the blade section having an area ab as shown in 
figure 1. Let the temperature of this section be TQ. The heat 
received in time dt through the two surfaces of unit span that are 
exposed to the gas stream is 


aq = (hg + hy)b(T& & p - Ty) at (1) 
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This heat dQ will increase the temperature of the blade strip accord- 
ing to the equation 


ат 
dQ = p bac | St. dt (2) 


Equating these two values of dQ (equations (1) and (2)) and solving 
for dt yield 


pac am 

= u: ШЫ» 22223 

dt E BAT тт (5) 
5 g,e,t m 


The boundary conditions for equation (3) are established when t = 0, | 
Tm == Tin, i = Tg,e,i end when t ze, Tm = Im,f ex Ig,e,f' 
Рша Ст 


If the term ұ---ұ- of equation (5) 1s considered constant and is 


called т, equation (3) can be integrated to yield 


Tg,e,f ` Tm -t/T 
m тр = e (4) 
g,e,t m,i 


The symbol T has the dimensions of time and is usually known as the 
time constant. The time constant is an indication of the rate at which 
the blade temperature will respond to changes in the combustion-gas 
temperature. If the time t in equation (4) is equal to the time con- 
stant T, the ratio (Tg e,f - т) Ты, e; > Tn, i) will equal approxi- 
mately 0.37. This means that the time constant represents the time 
required for the blade temperature to increase or to decrease 63 percent 
of the over-all difference between the initial and final blade tempera- 
tures. Equation (4) permits the calculation of the blade-metel temper- 
ature Ty, at any time t when the initial and final effective gas tem- 
perature, the blade-met&l density and specific heat, the local gas-to- 
blade heat-transfer coefficients, and the blade geometry are known. 


Case 2 - finite biade-metal thermal conductivity. - The more exact 
cage where the blade-metal thermal conductivity is finite is now con- 
sidered, so that the magnitude of the error that is introduced when 
blade temperatures are calculated with equation (4) can be evaluated. 
For the analysis, as stated previously, a strip of area of a blade sec- 
tion of unit depth was so chosen that no spanwise or chordwise heat flow 
would occur at the location. Consider such a strip of area ab of the 
turbine-blade section as shown in figure 1. A heat balance is made on 
the element of length dx of the strip of area ab. The plane O-O is 
so located that there will be no heat flow across this plane. The 
distance from the plane О-О to the suction surface is 7. The dis- 
tance from the plane 0-О to the element dx is x. The rate at which 


NACA RM Е5ЗА19 


heat is transferred from the element dx Беа the boundary &t а 
distance x from the plane О-О is 


am, 
kno Sx | (5) 


Lhe rate at which heat is transferred through the boundary at a distance 
x + dx is 


КЕБЕ ы " 


The net rate at which heat is gained by the element is 


Өт, 
стото 55” ( 7) 


When this new rate of heat gain is placed equal to the difference 
between the rate of gain and the rate of loss from the element dx the 


result is 


= x Ta 


шр ge = mu cce (8) 
Equa- 
tion (8) can therefore be written as барш 
M a^. Ж 
St 7 9m See 
oT 


Since the origin plane О-О was chosen as a plane where s= = O, 


the plane 0-0 would fall at the center of the strip only if the heat- 

transfer coefficients on opposite surfaces of the blade were equal. For 
a turbine blade these heat-transfer coefficients are not equal, and the 

plane О-О is located near the blade surface with the lower heat- 


transfer coefficient. For this reason, two solutions of equation (9) are 


needed, one for the region between О-О and the suction surface of the 
turbine blade and the other between the plane О-О and the pressure 
surface. 


Consider first the region between О-О and the suction surface. 
In order to obtain a solution of equation (9), it is convenient to 
introduce а new variable @=T - T,. in terms of this variable 0, 
g,e, f 
equation (9) becomes 
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The boundary conditions for equation (10) are 


(1) when x= O, S0/Əx = 0 


(2) when х= 1, - hg, = ky (38). 


(5) when о, Ө->0 


Assume that @ can be expressed as the product of two variables X 
and Y and that X is a function of x alone and Y is a function of 
t alone. The substitution of these valves in equation (10) and the 
solution of the equation as presented in appendix C yleld 


9 -Мкёа ё 
Ө = Се сов АКх (11) 


Asl 


From this equation, it can be seen that the value of the time constant T 
for the case of finite metal conductivity is defined by 1/A°K@a,. The 
integration constants C and K of equation (11) must be chosen to 
satisfy the three boundary conditions and also to make the two strip 
sections from the plane О-О to each surface of the blade compatible. 
To check boundary condition (1) differentiate equation (11) so that 


ш 2 
oF ac 
3% = < -KACe sin AKx (12) 


Thus, boundery condition (1) is satisfied for all values of C and K. 
Boundary condition (3) is also satisfied ot values of C and K, 
-N Keo, 


inasmuch as the exponential factor e approaches zero as % 
becomes infinite. Substitution of equation (11) in boundary condition 
(2) yields 


Bu E^ 
- hg < e Cg cos AKgl = “km L 1 e ACgKg sin AKsl (13) 


Equating corresponding terms in equation (13) results in 


hal 
— = AKgl ten AK,2 (14) 


Em 
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A transcendental equation of the same type as equation (14) is discussed 
in reference 4, where it is shown that there are an infinite number of 
roots. Reference 5 contains a table of the first six roots of this 
equation. For the case of the turbine blade being considered here, only 
the first root will be considered. Calculations using additional roots 
showed & negligible effect on the solution. Thus Kg is determined by 
the ratio bet /Ky: For the section of the strip between the plane 0-0 


and the suction surface, Og will be given by the relation 


Keen 


where K, is defined by equation (14). The section of the strip between 
the plane О-О and the pressure surface can be treated in а similar 
manner. өр will be given by 


On = Cse сов Kax (15) 


өр = o e Tont сов Kx (16) 


and. Ко is defined by the relation 
a-1 
a = K. (a-1) tan Ky(a-1) (17) 


If these equations are valid, they should give the same value for @ at 
the plane 0-О at all times. That is, 


2 
„Кол -Kamt 
баты та Ут (18) 


To ensure this equality at all times, C4 must equal Cp and Kg must 
equal Кр. The value of C is determined from equation (11) when t 
and x equal zero. The value of K is determined by solving equa- 
tions (14) and (17) simultaneously for К and 7. 


Coated Blade 


When a low thermal conductivity coating is applied to the turbine 
blade, its effect appears in equation (4) or (11) as a change in the 
blade-surface heat-transfer coefficient. In equation (4) the heat- 
transfer coefficient appears in the value of т. In equation (11) it 
appears in the value of K. The expression for the heat-transfer coef- 
ficient that should be used when an oxide coating is applied is pre- 
sented in reference 6. This equation, in tbe notation of this report, 
is 


1 1 c c 
m = h T = (19) 
Aa et M 
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It can be shown from the definition of the time constant and equa- 

tion (19) that the change in the time constant when a coating is applied 
is proportional to both бе/ Ке and h. Since the change in time con- 
stant is proportional to h, the use of low thermal conductivity coat- 
ings has a greater effect in applications where high heat-transfer coef- 
ficients exist, such as in rocket motors (ref. 1). 


APPARATUS AND INSTRUMENTATION 
Test Engine 


The experimental portion of the investigation was conducted on & 
production turbojet engine which had a centrifugal compressor and a 
single-stage turbine. The turbine rotor blades were standard service 
blades made of S-816 alloy. 


Fuel and Fuel Additives 


The fuel used for this investigation was MIL-F-5624A (grade JP-4). 
Two mixtures of fuel and fuel additives were used; one contained 6 per- 
cent tributyl borate by weight and the other contained 1.2 percent 
silicone oil, SF-99, by weight. These mixtures produced & l-percent 
concentration of boric oxide or silicon dioxide, respectively, in the 
produets of combustion. These additives were chosen for the investiga- 
tion because of their availability in the quantities desired. 


Instrumentation 


The engine thrust was measured with & calibrated air-celi thrust- 
meter. The engine speed was measured with a chronometric tachometer. 
The turbine rotor and stator blades were instrumented with 30-gage 
chromel-alumel thermocouple wire at the locations shown in figure 2. 
The blade thermocouples were made by threading the thermocouple leads 
through two-hole ceramic tubing. The ceramic tubing was inserted in a 
1/16-inch-outside-dismeter Inconel tube. The end of the Inconel tube. 
was spun over to seal the thermocouple into the tube. These thermo- 
couples were installed in radial grooves cut into the blade surface. 
When the thermocouples were located in the blade, the grooves were filled 
with a nickel-chromium brazing alloy which was then smoothed off and 
faired with the blade surface. 


The thermocouples shown in figure 2(a) were located in three dif- 
ferent blades so that no single biade would be excessively weakened by 
the grooves that were necessary for the thermocouple installation. 
Thermocouples 1 and 2 were located in one blade, 3 and 4 were located 
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in & second blade, and the thermocouples giving the spanwise tempera- 
ture distribution at the blade leading edge were located in a third 
blade. The leads from these thermocouples were fastened to the rear 
face of the rotor and then led through the drilled turbine and compres- 
sor shaft to the front of the engine where they were connected to a 
slip-ring thermocouple pickup. A similar type thermocouple system is 
described in greater detail in reference 7. Two of the stator blades 
were instrumented as shown in figure 2(b) to obtain steady-state meas- 
urements of the spanwise temperature distributions at the leading edge 
of the blades and to obtain stator-biade transient temperature data. 
The stator-blade temperature distributions, which were obtained only 
during steady-state conditions, are an indication of the radial 
combustion-gas temperature distributions leaving the two burner liners 
upstream of the two instrumented stator blades. The steady-state tem- 
perature data were obtained with an indicating potentiometer, and the 
transient blade temperature data were obtained with recording potenti- 
ometers that drew graphs of temperature against time. Samples of these 
graphs are shown in figure 5. 


EXPERIMENTAL PROCEDURE 


The transient blade temperature data used in this report were 
obtained at static sea-level test conditions. The data were obtained by 
accelerating the engine from 8000 to 11,500 rpm (rated engine speed) and 
recording the turbine rotor-blade temperatures during the transition 
between these two operating conditions and until the blade temperatures 
reached equilibrium at the new conditions. The accelerations were made 
from the initial engine speed of 8000 rpm to the final engine speed of 
11,500 rpm, because,for the turbojet engine used, acceleration through 
this speed range resulted in the largest difference in values of initial 
and final combustion-gas temperatures without exceeding rated engine 
speed. The engine was not accelerated above rated engine speed, because 
it was believed that coating thicknesses would not be sufficient to 
provide a blade temperature lag great enough to avoid blade failure at 
engine speeds above the rated value. Deceleration as well ав accelera- 
tion data were recorded in this investigation, so that a comparison 
between the theory presented in this report and the experimental results 
could be made for two different types of transient heat flow. Several 
engine accelerations and decelerations were made with each fuel additive 
mixture and also with the JP-4 fuel to obtain check data and also because 
there were not enough recording potentiometers available to record all 
blade thermocouple temperatures at once. It was therefore necessary to 
repeat test runs with different combinations of thermocouples connected 
to the potentiometers. 


The engine was operated first with no additive in the fuel, and data 
were obtained for eleven engine accelerations and seven decelerations to 
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obtain a reference for comparisons with the coated-blade data. During 
these runs time-temperature relations were obtained for thermocouple 
locations 1, 2, 3, and 4 on the rotor blade and for location 9 on the 
stator blade. Silicone oil was then added to the fuel, and the engine 
was operated at rated engine speed for approximately 1/2 hour to provide 
for the coating build-up before data were taken. Eight accelerations and 
eight decelerations were made to obtain transient blade temperature data 
for the same thermocouple locations that were considered during the 
uncoated-blade operation. After this period of running, the engine was 
disassembled and the engine parts were inspected for evidences of 
deposits of silicon dioxide. Measurements of the coating thickness were 
made where thickness and uniformity of the coating permitted. The coat- 
ing was then cleaned from the engine parts and the engine was reassen- 
bled. Six accelerations and six decelerations were made with the 
cleaned engine with JP-4 fuel used to check the first transient uncoated- 
blade temperature data that were obtained. The tributyl borate was 

then added to the fuel and the engine was operated for approximately 

1/2 hour at rated speed to provide for the coating build-up. Eight 
accelerations and eight decelerations were then made to determine the 
effect of this coating on the blade temperature-time relation. Follow- 
ing this running, the engine parts were again inspected to determine the 
thicknesses and the locations of the coatings. 


Bach time the engine was started and before time-temperature data 
were obtained, the engine was operated at an engine speed of 11,500 rpm 
and the turbine-blade temperature was set at 1365° F by adjusting the 
tail-pipe nozzle opening. This procedure produced transient blade tem- 
perature data that were always between the same temperature limits and 
also gave about the same engine combustion-gas weight flow and thus 
about the same blade-surface heat-transfer coefficients. Before addi- 
tives were mixed with the fuel, the engine thrust was measured to obtain 
a base value. During the fuel-additive investigations the engine thrust 
was measured to detect any change in the engine performance caused by the 
coatings on the turbine blades and other engine parts. 


CALCULATION PROCEDURE 
Theoretical Determination of Blade Time Constants 


The magnitude of the time constant, as indicated in the analysis, 
determines the rapidity with which the blade-metal temperature responds 
to changes in combustion-gas temperature. Therefore, to determine the 
effectiveness of a coating, the terms which influence the magnitude of T 
must be evaluated. For an uncoated blade, these terms are the gas-to- 
blade heat-transfer coefficients, blade geometry, and blade-metal con- 
ductivity, density, and specific heat. In addition, two other terms, 
coating conductivity and coating thickness, must be considered when the 
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blade is coated. Іп order to evaluate these terms, combustion-gas-flow 
conditions of a typical turbojet engine (ref. 8) were assumed. Also, 
in order to simplify the calculations, only a section at the blade 
midchord was considered. 


The procedure followed in determining the terms of the time con- 
stant was as follows: The velocity distributions along the surface of 
the blade profile at the midspan were obtained by the use of the stream- 
filament theory described in reference 9. The local velocities at the 
midchord surfaces of the blade with the conbustion-gas properties eval- 
uated at final blade-metal temperatures were then used to determine the 
local convective heat-transfer coefficients on the suction and pressure 
surfaces by methods presented in reference 10. The values of heat- 
transfer coefficients thus calculated are for uncoated blades. 


The effects of the coating were included by evaluating the over-all 
heat-transfer coefficients with equation (19) of this report. Since the 
coatings applied to the blade would probably be thin relative to the 
blade thickness (ref. 3), the effect of the coating on the velocity 
distribution around the blade periphery will in most cases be negligible, 
and it is so considered in the calculations of the convective heat- 
transfer coefficients in this report. However, for cases where the 
coating is thick, the blade profile may be altered enough to change the 
velocity distribution around the blade, with the result that the con- 
vective heat-transfer coefficient h would be different from that 
obtained for the uncoated blade. It would therefore be necessary to 
determine the velocity distribution and the heat-transfer coefficients 
for the modified profile before equation (19) could be used to determine 
the value of the over-all heat-transfer coefficient hyp. 


The magnitudes of T were calculated for both engine accelera- 
tions and engine decelerations by evaluating the metal density and 
specific heat at the final metal temperature. The blade-metal thermal 
conductivity was elther considered infinite or evaluated at the final 
metal temperature. 


Theoretical Determination of Effectiveness of Coatings 


The theoretical effectiveness of coatings was determined by calcu- 
lating, from equation (4), the blade-metal temperature variation with | 
time for various combinations of coating thickness and coating conduc- 
tivity necessary to produce a desired lag in blade-metal temperature. 
The value of the time constants used in equation (4) was calculated by 
assuming the blade-metal conductivity infinite and the same coating 
thickness on each surface of the blade, by neglecting the effect of 
coating thickness on the convective heat-transfer coefficients, and by 
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using values of over-all heat-transfer coefficients determined by equa- 
tion (19) for various ratios of coating thickness to coating 
conductivity. 


Method of Determining Blade Time Constants from Experimental Data 


Determining the effectiveness of the coatings obtained in the 
experimental investigation required accurate values of the time constant 
because of the relatively thin coatings deposited on the blade. The 
calculation procedure, therefore, is so arranged as to deduce the best 
possible value of T from the time-temperature data. Typical traces of 
the time-temperature variation at the midchord of a rotor blade as 
obtained with recording potentiometers are shown in figure 3. 


By assuming the simple case of infinite blade-metal thermal conduc- 
tivity (which, as will be shown later, varies only slightly from the 
more complex case of finite blade-metal conductivity) and the engine 
conditions as changed instantaneously, the theoretical equation of a 
blade temperature response with time was given in the analysis section 
as 


T - Т -t/T 
g,e,f n ue / (4) 
Ig,e,f - Tm,i 
When the case of an uncoated blade is considered in which the bilade- 
metal conductivity is infinite, the value of T is given for the 
midchord region of the blade as 
р. ас 
T = P m_ 
hs + bp 


where, with the assumption that the change from one temperature level 
and one set of flow conditions to another is instantaneous, values of 

the convective heat-transfer coefficients hg and hp are fixed and Т 
is a constant. However, in the experimental investigation, approximately 
6 to 15 seconds were required to change engine conditions, with the 
result that the heat-transfer coefficients varied and thus the value of 

T also varied during the first 6 to 15 seconds. Consequently, the value 
Of T ів determined from the data obtained after final combustion-gas- 
flow conditions are established or approximately 10 seconds after the 
time the engine speed is changed. 


If equation (4) of the analysis is written in the logarithmic form, 


In(Tg,e,f - Tm) = 1n(Tg,e,f - Тш,1) - t/t (20) 
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it can be seen that, since In(Tg,e,r - Tm,i) is а constant, the data 


can be represented on a semilog plot as a straight line, if the temper- 
ature difference (Ty e,f - Т) is plotted on a logarithmic scale and 


the time on an equal-parts scale. Equation (20) also indicates that the 
negative reciprocal of the slope of the line on this plot represents the 
value of the time constant T. The data obtained while the engine con- 
ditions axe changing will not fall on a straight line, because of the 
variation of the blade-surface heat-transfer coefficient during this 
period; however, after the engine combustion-gas-flow conditions are 
established, the data should fall on a straight line, as is illustrated 
by & typical plot of the dat& shown in figure 4. The intercept of the 
straight line on this type plot &t zero time gives & fictitious value of 
the temperature difference Tg,e,f - Tm,i that would satisfy the theo- 
retical conditions of instantaneous change in engine conditions. 


The value of T found by this graphical method depends to some 
extent on judgment in fairing a line through the data and in selecting 
a value of Tg,e,f- The first value of Tg,e,f usually tried is that 
which appears to be the final steady reading from the time-temperature 
traces (fig. 3). However, for some runs, after apparently reaching & 
final value, the final metal temperatures indicated by the traces 
drifted a few degrees so that the value obtained made the data incom- 
patible with theory. A small constant error in the value of Tg,e,f, 
when included in а plot such as figure 4, makes the data deviate from 
the straight line indicated by theory. The deviation from the straight 
line becomes greater when time increases and Tg,e,f - Tm десгеавев. 
The constant error in Tg,e,f then becomes an increasingly large per- 


centage of the temperature difference Tg,e,f - Tm. Because of this 
known inaccuracy in the measurement of Tg,e,f> the observed value was 
adjusted by a constant for any one run (never more than the instrumen- 
tation error) to produce a straight line as predicted by theory when the 
data is plotted as shown in figure 4. 


Tn order to obtain as accurate a time constant from the data as 
possible, a method of least squares presented in appendix D was used to 


refine the approximate graphical results. The graphical results, because 


of their dependence on the observer's judgment in selecting appropriate 
values of Tg,e,f and in fairing a line through the data, for some runs 


indicated a difference of several seconds in the value of T as compared 


with values obtained with the least-square method. 
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RESULTS AND DISCUSSION 
Experimental Time Constants 


Experimental time constants were determined for each of the ther- 
mocouple locations on the rotor blade and for the leading edge of the 
stator blade with the method of least squares outlined in appendix D. 
The time constants were obtained by using the transient Бізде temperature 
даља for the time after the combustion-gas flow was established at the 
final conditions. The time constants obteined in this investigation for 
both accelerations and decelerations are summarized in table I. 

Table I(a) contains the time constants for the engine accelerations for 
the coated and the uncoated blades. There is some scatter of the time 
constants obtained from individual runs for each thermocouple location, 
and for this reason these values are averaged to obtain values for com- 
parison between the coated and the uncoated blades. An inspection of 
the average time constants for the rotor blades reveals that for all the 
thermocouple locations on the turbine rotor blades the coatings of both 
silicon dioxide and boric oxide formed by the fuel additives (silicone 
oil, SF-99, and tributyl borate, respectively) had no measureable effect 
on the temperature response of the turbine rotor blades. The variation 
between the average time constants for coated and uncoated blades is of 
the same magnitude as the variation obtained for individual runs made on 
the uncoated blades. This variation between coated and uncoated blades 
is therefore within the experimental accuracy of the apparatus, and the 
effect of the coatings on the rotor blades is negligible. 


The average experimental time constant for the leading edge of the 
uncoated stator blade was 21 seconds. The silicon dioxide coating on 
the stator blades increased this value to 29 seconds. This difference 
in time constant would result in about a 25° F decrease in the stator- 
blade temperature relative to the uncoated blade 60 seconds after the 
engine operating conditions were changed. A more complete discussion 
of the relation between time constants and the resulting blade tempera- 
ture is presented later in this discussion. The boric oxide coating on 
the stator blade resulted in a time constant of Z5 seconds, which was 
slightly higher than the uncoated-blade time constant of Zl seconds; but 
the difference between this value and that of the uncoated blade is 
within the accuracy of the experimental data. 


Table I(b) contains the values of the time constants obtained when 
the engine was decelerated from 11,500 to 8000 rpm.. These values of 
time constant show the same effects of the coatings as were shown in the 
acceleration data. The greatest effect was at the leading edge of the 
stator blade when the blade was coated with silicon dioxide. The decel- 
eration time constants were larger than those obtained for engine 
acceleration, because the time constant is an inverse function of the 
local gas-to-blade-surface heat-transfer coefficient which exists at the 
final condition (11,500 rpm for acceleration and 8000 rpm for 
deceleration). 


p -————— all 
16 ONF IDENTIAL NACA RM E53A19 


Characteristics of Coating Deposits on Blades 


Inspection of the turbine blades after using each of the fuel addi- 
tives showed that both coatings were very thin. This thinness was the 
reason for their ineffectiveness in decreasing the blade temperature- 
response rate. The silicon dioxide coating had a powdery appearance and 
Showed a tendency to flake off the blades. This coating was powdery 
because the engine combustion-gas temperature was less than the melt- 
ing point of silicon dioxide. The boric oxide coating was very thin 
and had & glassy appearance. The thickest coating was that of silicon 
dioxide on the stator blades. This inspection verified the results 
presented in table I, which shows that the leading edge of the stator 
blades had the greatest increase in the time constant. Measurements of 
the coating thickness were made for the siliton dioxide coating; but, 
because it flaked off easily and because the flow conditions over the 
blade were altered as the engine was shut down, the measured thicknesses 
are probably not the same as those that existed when the data were taken. 
Also, because the coatings were quite thin and powdery, it was difficult 
to obtain accurate measurements of the coating thickness. The boric 
oxide coating was so thin that it was impossible to determine the thick- 
ness that had been deposited on the rotor or the stator blades. The 
measured thicknesses of the silicon dioxide coating on the turbine rotor 
blades after the first hour of engine operation were: for the suction- 
surface leading edge, 0.0005 inch; midchord, 0.001 inch; and trailing 
edge, 0.005 inch. The pressure surface had а uniform 0.0005-inch coating 
over the entire surface. The leading edge of the stator blade had & 
0.015-inch coating. 


Comparison of Time-Temperature Curves Calculated from 
Experimental Time Constants for Coated and Uncoated Blades 


The average values of the experimental time constants for the stator- 
blade leading-edge thermocouple for the uncoated blade and the blade 
coated with silicon dioxide were used with equation (4) of the analysis 
to calculate curves of the blade temperature variation with time. The 
initial and final blade temperatures assumed for the calculation were 
1000° F and 1365° F, respectively. These two curves are presented in 
figure 5. The difference in the time constants for these two curves is 
8 seconds, and it is reflected in a blade temperature lag of about 25° F 
at a time 60 seconds after the operating conditions of the engine were 
suddenly changed. Similar time-temperature curves plotted for the tur- 
bine rotor-biade locations would show no difference between the coated- 
and the uncoated-blade curves, because of the small variations in the 
rotor-blade time constants which are presented in table I. 
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Comparisons of Theoretical and Experimental] Uncosted-Blade 
Time-Temperature Variations and Time Constants 


Two equations, (4) and (11), derived in the analysis of this report 
permit the calculation of the blade temperature at any time t after 
the engine operating conditions have been suddenly changed. Equa- 
tion (4) is based on the assumption that the thermal conductivity of the 
blade metal is infinite, while equation (11) is based on an assumption 
of a finite blade-metal thermal conductivity. These two equations were 
used to calculate the variation of the blade-metal temperature at the 
midchord, midthickness location of the uncoated turbine rotor blade for 
the conditions of this experimental investigation. No data from the 
experimental investigation were used in the calculation. The two cal- 
culated curves of the blade time-temperature relation are presented for 
engine acceleration and compared with an experimental curve in figure 6. 
The values used in obtaining the calculated curves are presented in the 
following table: 


Blade-metal properties Blade-surface heat- 
transfer coefficients, h, 
Btu/(sec) (sq ft) (°F) 


Thermal Suction- Pressure- 
conductivity, Ems surface surface 


Btu/(sec) (£t) (OF) midchord midchord 


556 0.062 0.015 


The surface heat-transfer coefficients were calculated and the blade- 
metal properties wére evaluated as described in the CALCULATION PROCEDURE 
of this report. The calculated value of the heat-transfer coefficient 
at the suction surface is much larger than that at the pressure surface, 
because, for the combustion-gas-flow conditions at an engine speed of 
11,500 rpm, the flow at the suction surface was turbulent while that at 
the pressure surface was laminar. The calculated time constants for the 
cases of infinite and finite blede-metal thermal conductivity were 20 
and Zl seconds, respectively. These time constants were used in equa- 
tions (4) and (11) to calculate the blade time-temperature relation. 

The time-temperature curves which were obtained with these two values of 
time constant are presented in figure 6. 


Considering the blade-metal thermal conductivity infinite greatly 
simplifies the calculation of the blade-metal temperature variation and 
as shown in figure 6 results in & maximum difference of 6° F between the 
curves calculated by the two equations (4) and (11). Similarly, calcu- 
lations of the magnitude of T made for the case of engine deceleration 
indicated small effect due to considering the conductivity to be either 
finite or infinite. The computed values of T for cases of infinite and 
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finite metal conductivity were 29 and 31 seconds, respectively. Тһе 
larger time constant was obtained for the case of a finite blade-metal 
thermal conductivity which results in a finite time for the heat to flow 
to the midthickness location. Since the time constants obtained for the 
cases of finite and infinite blade-metal thermal conductivity agree, the 
latter, being less complex, was used in the calculations of the tran- 
sient blade-metal temperatures presented in this report. The values 
would not agree for a blade that was considerably thicker than the blades 
used in this investigation or for blades of very low thermal conduc- 
tivity, such as ceramic blades. 


The calculated values of time constant at the midchord location and 
the experimental values obtained from the data after final engine 
combustion-gas-flow conditions were established agree very well. The 
calculated values as defined in equation (4) for the case of infinite 
thermal conductivity were 20 seconds for acceleration and 29 seconds for 
deceleration. The average experimental values were 20 seconds and 
$0 seconds for acceleration and deceleration, respectively. Therefore, 
& time-temperature curve calculated with the theoretical value of the 
time constant would be the same as one calculated with the experimental 
value. However, when the calculated time-temperature curve for an 
engine acceleration is compared with a plot of the blade temperature 
дада, against time, exact agreement between the calculated and the data 
curves cannot be obtained, because the engine operating conditions can- 
not be changed instantaneously from 8000 to 11,500 rpm as is assumed in 
the derivation of equation (4). A comparison between the curve cal- 
culated with equation (4) and the data curve plotted from run 14, 
table I(a), is shown in figure 6. The maximum difference between the | 
two curves is about 40° F. The time required to change the engine 
operating conditions from 8000 to 11,500 rpm for this run was 7 seconds. 
The magnitude of the disagreement between the calculated and the data 
curve is a function of the time required to change conditions. The 
shorter the time required to change the operating conditions, the better 
will be-the agreement between the calculated and the data curves. 


Determination of Desired Coating Properties 


The comparison between experimental and theoretical determination 
of turbine rotor-blàde time constants shows that equation (4) is reliable 
for calculating the uncoated-turbine-blade: time constant. It was not 
possible to check the reliability of equation (19) with the data of this 
investigation, because the coating thicknesses obtained could not be 
measured accurately. However, in reference 1 data are presented for the 
temperature reduction of rocket motor walls when they are coated with 
several different refractory materials. Coating-thickness data were 
presented only for zirconium dioxide coating. By using equation (19) 1% 
was possible to predict the reduction of heat transfer into the rocket 
motor walls within an accuracy of 15 percent. Therefore, equations (4) 
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and (19) will be used to determine the effect of different ratios of 
coating thickness to coating conductivity 6,/k, on the turbine-blade 


temperature response. A series of calculated time-temperature curves 
for a range of 6¢/ke (from О to 50 (sq ft)(sec)(°F)/Btu) is presented 
in figure 7. These curves have been calculated for the midchord loca- 
tion on the blade, with the assumption that the coating thicknesses on 
the suction and pressure surface are the same. 


If, for the purpose of illustrating the shortcomings of dynamically 
applied oxide coatings, an arbitrary lag of the coated-blade temperature 
of 100° F relative to the uncoated blade is desired 60 seconds from the 
time the engine operating conditions are changed from 8000 to 11,500 rpm, 
it can be seen from figure 7 that a value of 5./Е. equal to about 


57 (ва ft)(sec)(9?F)/Btu is required for the engine combustion-gas-flow 
conditions used in this investigation. The silicon dioxide coating 
thickness required to produce this lag in turbine rotor-blade tempera- 
ture could not be calculated, because the thermal conductivity of the 
powdered silicon dioxide deposited on the turbine blade in this inves- 
tigation was unknown. Available literature gives the thermal conduc- 
tivity only of fused silicon dioxide, which is higher than the powdery 
form that was obteined on the turbine blades. However, by using the 
measured values of the coating thicknesses obtained in this investi- 
gation with silicon dioxide, approximately 0.001 inch, it is possible 
to calculate the coating conductivity required to produce this arbi- 
trarily desired lag in metal temperature. A coating thermal conductivity 
of 0.008 Btu per hour per foot per °F would be needed for a coating 
0.001 inch thick. For a coating thickness of 0.010 inch, the maximum 
considered permissible from aerodynamic considerations of this turbine, 
& coating thermal conductivity of 0.08 Btu per hour per foot per Ор 
would be required to produce the temper&ture lag of 1009 F. А survey of 
literature on possible oxide coatings indicated that there were no mate- 
rials available with thermal conductivities as low as 0.008 Btu per hour 
per foot per OF but that materials do exist with thermal conductivities 
in the range of 0.08 Btu per hour per foot per °F. One such material 

is antimony oxide, which Тав a conductivity of 0.068 Btu per hour per 
foot per ОЕ. This oxide is a possibility for a blade-coating material 
if it can be dissolved in the engine fuel and if a method can be 
devised to cause the coating to build up to a thickness of 0.010 inch. 


Deposits Formed on Engine Parts 


Only the silicon dioxide deposits in the engine are presented in 
photographs, because the boric oxide deposits were so thin for the con- 
centration of the additive used as to be negligible. The photographs give 
a clear idea of the types of coating encountered and the problems of 
engine operation that may result from these coatings. The silicon 
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dioxide coating that formed on the turbine rotor blades is shown in 
figure 8. The uncoated blades are shown in figure 8(a), and the blades 
after one hour of engine operation in figure 8(b). After this first 
hour of operation, the coating covered the entire blade. The coating 
thicknesses were measured with the coating in this condition. After an 


additional 15 hours of operation, the turbine rotor blades were coated 


as shown in figure 8(c). The coating did not adhere to the blades as 
well as it did after the first period of operation, and, as shown in 
figure 8(c), most of the coating on the blades had chipped off. The 
chipping probably resulted from loosening of the coating when the blades 
cooled after the first engine inspection shut-down. When the engine 
was started again the loosened coating was washed off. 


The turbine stator blades are shown in figure 9; part (a) before 


engine operation and part (b) after 22 hours of operation. The deposit. 
on the stator blades was heavier than on the rotor blades, but because 
of its thickness the coating appeared to flake off after it reached a 
given thickness and then to build up again. The complete washing off of 
the coating appears to be a characteristic only of the rotor blades, for 
no such complete washing was evident on the stator blades. 


Photographs of the burner domes after 2$ hours of engine operation 


are shown in figure 10. The burner domes were coated with a thick, 
spongy coating of silicon dioxide and carbon. The thickness of the 
deposit was about 1/8 inch. There was also a build-up of deposit on 
the fuel nozzles which in time could conceivably alter the shape of the 
fuel spray and decrease the efficiency of combustion. The burner liners 
also had a deposit about 0.010 inch thick on the inside surface that was 
similar to that in the burner domes. | 


Only deposits of silicon dioxide in the engine have been discussed, 
because they were much greater than the boric oxide deposits formed by 
the additive tributyl borate. The boric oxide deposits were so slight 
as to be immeasurable. A characteristic of tributyl borate is its 
affinity for water. When it combines with water, boric oxide precipi- 
tates from the additive solution and thereby renders the additive inef- 
fective for forming coatings on the turbine blades. This precipitate 
can also cause engine operating difficulties by clogging fuel filters; 
and, if the precipitate forms while the fuel is in the engine, diffi- 


culties such as the plugging of fuel nozzles can occur. During the oper- 


ation of the engine in this investigation the fuel nozzles became 
Plugged with the boric. oxide precipitate which formed when the fuel 
remaining in the fuel nozzles combined with the moisture in the air 
after the engine was shut down for an inspection of the turbine blades. 
Photographs of three typical plugged fuel nozzles are compared with 
clean nozzles in figure 11. Nozzles number 4 and 7 were partially 
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plugged, and nozzle 11 was completely plugged. The plugging was severe 
enough to make it impossible to start the engine. After this difficulty 
was encountered, engine shut-downs were preceded by purging the fuel 
system of the engine with straight JP-4 fuel to clean the additive from 
the fuel nozzles. This procedure was effective in eliminating the fuel- 
nozzle plugging difficulties. 


Effect of Fuel Additives on Engine Performance 


The use of either SF-99 silicone oil or tributyl borate as an addi- 
tive in the fuel for a period of approximately 2 hours of engine oper- 
ation did not cause a measurable change in the engine thrust, as com- 
pared with an engine operating on JP-4 fuel with no additive. There 
also was no change in the burner-outlet spanwise temperature distribution 
resulting from the addition of the silicone oil, SF-99, or the tributyl 
borate to the fuel. The only change in the engine performance was an 
increase in the specific liquid (fuel plus additive) consumption. This 
change was a result of depression of the heating value of the JP-4 fuel 
by the low heating values of the additives, which made a higher fuel- 
air ratio necessary to obtain the required turbine-inlet temperature at 
rated engine speed of 11,500 rpm. 


SUMMARY OF RESULTS 


An analysis was made to permit the prediction of the temperature 
response of coated and uncoated uncooled turbine rotor blades after a 
sudden change in the combustion-gas temperature. An experimental inves- 
tigation was made on & turbojet engine to check the analysis and to 
determine the effectiveness of dynamically applied coatings of silicon 
dioxide and boric oxide. The fuel additives used to produce these coat- 
ings were silicone oil, SF-99, and tributyl borate, respectively. The 
mixtures of these additives with the fuel produced a l-perceni concen- 
tration of the coating material in the products of combustion. The 
results of this investigation are as follows: 


1. Neither the silicon dioxide nor the boric oxide coating on the 
turbine rotor blades was effective in retarding the rate of turbine- 
blade temperature response. 


2. The analysis permitted calculation of the experimental transient 
turbine-blade temperatures of this investigation with a maximum error of 
about 40° F. 


2. For the coating thickness obtained with silicon dioxide at the 
midchord location (0.001 in.) the analysis indicated that a coating 
thermal conductivity of 0.008 Btu per hour per foot per °F would be 
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needed to produce a blade temperature lag of 100? F 60 seconds &fter the 
engine operating conditions were changed from 8000 to 11,500 rpm. A 
survey of literature showed that no oxide coatings having conductivities 
this low are available. 


4. For the maximum coating thickness (0.010 in.) that can be per- 
mitted on the turbine blade from aerodynamic considerations, the anal- 
ysis indicated that, for the metal temperature lag of 100° F 60 seconds 
after the engine operating conditions were changed, a coating thermal 
conductivity of 0.08 Btu per hour per foot per ОЕ was needed. A survey 
of literature indicated that only few oxides have thermal conductivities 
this low. Antimony oxide is one that falls in this range and is a pos- 
sibility for a blade-coating material if it can be used with the dynamic 
coating technique and if it will deposit on the turbine-rotor blades to 
& thickness of 0.010 inch. 


5. Engine operation for a period of approximately 2 hours with each 
fuel additive and the resulting deposits on the engine parts had no 
measurable effect on the engine thrust. 


Lewis Flight Propulsion Laboratory 
National Advisory Committee for Aeronautics 
Cleveland, Ohio 
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SYMBOLS 


The following symbols are used in this report: 
area of blade leading-edge section, sq ft 
length of element at midchord of blade, ft 
width of element at midchord of blade, ft 
integration constant 

specific heat, Btu/(1b)(°F) 


local gas-to-blade-surface heat-transfer coefficient, 
Btu/ (sec)(sq ғ+) (CF) 


average g&s-to-blade heat-transfer coefficient over length L of 
blade leading edge, Btu/(sec)(sq ft)(9F) 


integration constant 

thermal conductivity, Btu/ (sec) (ft) (°F) ; 
length of leading-edge periphery exposed to gas stream, ft 
distance from plane O-O to suction surface, ft 


intersection of suction surface and plane which separates leading 
edge from remainder of blade 


intersection of pressure surface and plane which separates leading 
edge from remainder of blade 


direction normal to boundary MN 


plane through blade element where ат /àx = 0 


heat flow, Btu 
length of leading-edge boundary, ft 
temperature, OF 


time, sec ` 


=» — h.. < ш ы F ж 
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X function of x 

x distance from plane O-O to incremental element of blade midchord 

section, ft | 

funetion of t 

о thermal diffusivity k/pc, sq ft/sec 

A correction to Та,е,Р,о 

5 thickness, ft 

€ correction to Tg 

e (Tee, А», OF 

À index of summation 

р density, lb/cu ft 

T time constant pac/h, вес 

Subscripts: 

c coated 

e effective 

f final 

g combustion gas 

1. initial 

m metal 

O values obtained by approximate graphical method 

Tp pressure surface 

S suction surface 

E applied to integration constants 
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APPENDIX B 


ANALYSIS OF TRANSIENT METAL TEMPERATURES AT LEADING OR 
TRAILING EDGES OF TURBINE BLADES 


The equations presented in this section are derived only for the 
leading edge of the turbine blade, but they also may be applied to the 
trailing edge. The only differences between these two locations are the 
blade geometry and the gas-to-blađe heat-transfer coefficients. 


Uncoated Blade 


Consider a section of the leading edge of a turbine blade as shown 
in figure i. This section has unit span, area A, and length L, of 
the periphery exposed to the exhaust gas. The plane M-N separates the 
leading-edge section from the rest of the blade. It is chosen such that 
the chordwise heat flow from the leading edge is perpendicular to the 
plane at all points. The direction of heat flow is n, and an element 
of Бірде in the direction of heat flow is dn. 


The net heat flow into the area A through the boundaries must equal 
the heat stored in the area. Such a heat balance results in the equation 


L=L 
Мае, ғ 7 Tm)GL +k, КЕ; ds = TUR E uy (81) 
from which 

а = Pm n Tm (B2) 

I=L N 

WT, е,ғ - T,,)aL + «| 

І=0 М 
= (83) 


(Tg,e,f - Tm) han + PIT 
І=0 
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If the boundery M-N is chosen so that (am, / dn)ds is zero or negligible, 


N 
ат | 
the term K, | (ша) ds can be neglected (the case of infinite blade- 
M 


I=L 
metal thermal conductivity). Also, if h(dL) can be evaluated 
L=0 
in terms of an average heat-transfer coefficient h for the length L 
of the leading-edge surface exposed to the gas stream, then 
L=L 
h(àL) = hL, and equation (B3) can be written 
L=0 


ыз РЕСЕ ®lm o (ва) 
(Т.е, - Ty) 


dt = 


The time constant for the leading edge of the blade then is defined as 


T = Pot m (B5) 
hI 


and integration of equation (B4) will again produce equation (4) 
Tg,e,f - Tm іш е-5/т (4) 
Tg,e,f 7 Tm,i 


If the heat flow across the boundary M-N cannot be neglected (the 
case of a finite blade-metal thermal conductivity), then the term 


N 
ат 
Km | (Ша) ds must be included in the time constant. The time con- 
M 
stant 1s then 


Pm Cm (вв) 


Тл [a s АНЕ ЕЕЕ 
hL(Tg,e,f - Tm) 


Coated Blade 


The effect of an oxide coating on the leading and the trailing edge 
of a turbine blade is handled the same as the coating on the midchord 
location. Equation (19) is used to calculate the over-all heat-transfer 
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coefficient for the leading edge of the coated blade. This coefficient 
is then used in equation (B5) or (B6) to determine the time constant, 

which is then used with equation (4) to determine the leading-edge blade 
temperature at any time t after engine operating conditions have been 


changed. 
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SOLUTION OF DIFFERENTIAL EQUATION FOR TRANSIENT BLADE TEMPERATURE WHEN 
BLADE-METAL THERMAL CONDUCTIVITY IS FINITE 
When the turbine-blade-metal thermal conductivity is considered 
finite in the calculations of the transient turbine-blade temperature 
after the turbine operating conditions have been changed, it is nec- 
essary to solve the differential equation 


36 ae 
557 “m 5-2 (10) 


х 


The solution of this differential equation is presented here. The 
boundary conditions for the equation are 


when x= 0, 06/dx = 0 
when x= 2, - hg@, = x, (32) 
2 s “8 | Ox/, 


when +оо, 0—0 


Assume that Ө can be expressed as the product of two variables X and 


XY, where X is а function of x and Y is a function of t. Thus 
Ө = X'Y (C1) 


Substitution of this value of Ө in equation (10) yields 


a 
ay d X 
at Х = Ont 21) (C2) 


or 


зар — (C3) 


Since the left term of equation (C3) is a function of t, and the right 
term is a function of x, both terms are constant and equal. This con- 


stant value can be written as -K^a,. Equation (C2) can now be written 
as two ordinary differential equations 


T TES (c4) 
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and 


2 
2 EA K (C5) 


These two equations can be rewritten &s 


T = -Kait (c6) 
and 

аёх 

"m K°X = O (ст) 


Equation (C6) can be integrated to produce 


in Y = Kat + in Y4 


or 


-Кёо, t 


Y = Y;e (C8) 


Equation (C7) is a simple harmonic equation and its solution can be 
written in а number of forms. The one best adapted to the present 
boundary conditions is 


X = X, cos Kx (C9) 


A solution of equation (Cl) can then be written as 


Ө = Y,X,e "m cos Kx (C10) 


The product X,Y, is a constant of integration to be determined from 


the boundary conditions. It is numerically equal to the temperature 
difference Ө at the plane О-О when t х= О and is defined ав а con- 
stant C. Because equation (C10) is only one solution of the linear 
equation and any combination of solutions of the equation is also s, 
solution, equation (C10) in its most general form can be written as 


- t 
0 S os "m cos AKx (11) 
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APPENDIX D 


METHOD OF LEAST SQUARES AS APPLIED TO APPROXIMATE GRAPHICAL RESULTS 
FOR DETERMINATION OF TIME CONSTANTS 


The procedure followed to obtain the time constants from the exper- 
imental data is presented herein. Equation (4) is solved for То and 
becomes | 


t/T 


T, = T - (T (D1) 


ш = Tg,e,f g,e,f ~ Tm,1)e 


Ihe quantities to be determined are Tg,e,p and т. These terms are to 


be given magnitudes which make the sums of the squares of the residuals. 
a minimum. А residual is the temperature difference (Tm,o - Tm), where 


-t/T 
Tmo = Tg,e,f,o - (Tg,e,f,0 - Тш,1)е % (D2) 


Because equation (D1) is not linear in Т, but transcendental, it will . 


be linearized by means of Taylor's theorem for two independent variables. 


Let 


Ta ЕЕ 


ш O 


Tg,e,f B Tg,e,f,0 + À 


where т. and Tee, f,o are the values found by the approximate graph- 
ical method, and € and A are the corrections to be made so as to 
minimize the sum of the squares of the residuals. Then by Taylor's 
theorem, 


or or 
Tn = Im о Қасты! + € e к... (D3) 
2 g,e,t O T O 


This equation is linear in A and є if these corrections are small 
enough that the terms involving the higher powers of them, such as 4f, 
€^, є2 , апа so forth, сап be neglected. А11 the terms in equation (р5) 
can be obtained from the data except the two unknowns A and €. 


The values of A and ç that render the sum of the squares of the 
residuals a minimum are given in reference ll as 


2785 


9822 
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Уве 0), 80) еә (ка) 
у) Ge) lal) 0 ДЕНЕ 


(D4) 


In equations (D4) the value of the blade-metal temperature Tm is 
obtained from the time-temperature data, while other terms are obtained 
by solving or differentiating equation (D2). By inserting the values of 
these terms obtained for each data point, the products and accumulations 
required in equation (D4) can be computed and the equations solved for A 
and €. The values of Tz, е,?,о and Tg obtained by the approximte 
graphical method then can be corrected by A and €, respectively, to 
obtain the corrected values of Tg,e,p and T. 
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TABLE I. - TABULATION OF BLADE TIME CONSTANTS 


| 


(a) Data obtained by accelerating engine from 8000 to 11,500 rpm. 


Time constant, T, sec 


Rotor-blade thermocouple location Stator-blade thermo- 
couple location 
Leading | Trailing Leading edge 


edge edge 


C i> WO [7 


ОО 0-30 


MIL-F -5624A (grade JP-4) fuel plus 1.2 percent by weight 
of silicone oil, БЕ-99 


МТІ.-Е-5624А (grade JP-4) fuel plus 6 percent by weight 
of tributyl borate 


Average 


denm 


55 


34 
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— 


TABLE I. - TABULATION OF BLADE TIME CONSTANTS - Concluded , 


(b) Nata obtained by decelerating engine from 11,500 to 8000 rpm. 


Time constant, T, sec 
Rotor-blade thermocouple location Stator-blade thermo- 
couple location 
Leading | Trailing |Midchord |Midchord Leading edge 
edge edge suction |midthickness 
GG C.a 


277 


-Е-5624А (grade JP-4) fuel 


оо чо лю Q EO HJ 


MIL-F-5624A (grade JP-4) fuel plus 1.2 percent by weight 
of silicone oil, SF-99 


MIL-F-5624A (grade JP-4) fuel plus 6 percent by weight 


of tributyl borate 


ЕЕЕ: 


== 29 Бс == 56 
Average |31 — | 29 | a | 32  — | 32 = 


2785 


55 
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Figure 1. - Turbine-blade areas considered in analytical development 
of blade time constants. 
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Cross section at 25 in. fram tip Cross section at 2 in. from tip 


а Thermocouple 


“RACAL” 


(a) Rotor blade. (b) Stator blade. 


Figure 2. - Locations of thermocouples on turbine rotor and stator blades. 
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Time, 10 вес per division———*- 


QU IUENTIAT, | 


(а) Acceleration; run 14. Thermocouple 
located at midchord, midthiekness of 
biade. 
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(p) Deceleration; run 9. Thermocouple 
ocated at midchord, midthickness of 
blade. 


Figure 3. - Typical time-temperature traces of 
uncoated-rotor-blade temperature obtained with 
recording potentiometers. 
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° ини located 
at midchord, mid- 
thickness of blade 


КП 


Temperature difference, Tg & p - Ty OF 


Time, t, sec 


Figure 4. - Typical plot of blade temperature data for determination 
of time constants. Acceleration; run 14. 
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Blade-metal temperature, Ta, OF 


1400 


1200 


Coated with silicon 
dioxide (based on 
average T = 29 sec, 
table I(a)) 

— — — — Uncoated (based on 

average T = Zl sec, 

table I(a)) 


1000 
0 


Time, t, вес 


Figure 5. - Comparison of calculated stator-blade time-temperature 
relations for uncoated blade and blade coated with silicon 
dioxide. Curves calculated with experimental time constants in 
equation (4); blade-metal initial temperature, 10000 F; blade- 
metal final temperature, 13659 F. 
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'd 1100 ? — — — — Theoretical curve; Т = Zl sec 
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" conductivity ) 
E —-—Lheoretical.curve; T = 20 sec 
= (considering infinite blade-metal 
x: conductivity) 
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О 20 40 60 80 100 
. Time, t, вес 


Figure 6. - Comparison of experimental and theoretical time- 
temperature relations at midchord, midthickness of uncoated 
rotor blade. Blade-metal initial temperature, 9729 F; 
blade-metal final temperature, 15650 F. 
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1400 


Ratio of coating thickness 
to conductivity, 8 /k s 


(sq ft)(sec)(9F)/Btu 


1500 
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u/7 
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Time, t, sec 


Blade-metal temperature, Т, OF 


Figure 7. - Theoretical variation of time-temperature relation at 
blade midchord with ratios of coating thickness to coating con- 
ductivity.  Blade-metal thermal conductivity considered infinite; 
blade-metal initial temperature, 10009 F; blade-metal final 


temperature, 15650 F. 
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С. 31586 
(c) Blades after 22 hours of operation. 


Figure 8. - Turbine rotor blades before and after operation with silicone oil, SF-99, added 
to the fuel. 
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(b) Blades after 21 hours of operation. c. 31587 > 
Figure 9. - Stator blades before and after operation with silicone oil, 5-99, added to the fuel. 


dat 


Figure 10. - Silicon dioxide and carbon deposits on burner domes after approximately 25 25 hours of operation with allicone 
oil, SF-99, added to the fuel. 
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Figure 11. - Fuel nozzles after approxima 
to the fuel compared with clean nozzles 
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